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Abstract 

The exact relation between the response function R(t, t') and the two time 
correlation function C(t,t') is derived analytically in the one dimensional ki- 
netic Ising model subjected to a temperature quench. The fluctuation dissi- 
pation ratio X(t,t') is found to depend on time through C(t,t') in the time 
region where scaling C(t,t') = f(t/t') holds. The crossover from the non- 
trivial form X(C(t,t')) to X(t,t') = 1 takes place as the waiting time t w is 
increased from below to above the equilibration time t eq . 
PACS: 05.70.Ln, 64.75. +g 
e-mail: lippiello@sa.infn.it, zannetti@sa.infn.it 

I. INTRODUCTION 

The time evolution of a system relaxating to equilibrium is characterized by two distinct 
time regimes: the off equilibrium transient for t < t eq and the stationary equilibrium evolu- 
tion for t > t eq , where t eq is the equilibration time. Normally, t eq is smaller than experimental 
times and one can actually observe equilibration. However, the existence of situations in 
which the opposite is true, namely where t eq exceeds by far any practical observation time 
as in glassy systems at low temperature, has contributed to arise a lot of interest in the 
off equilibrium relaxation regime. It turns out that quite a useful tool in the investigation 
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of these slowly relaxating systems, with and without disorder, is the fluctuation dissipation 
relation. 

For definiteness, let us consider a system in equilibrium at some temperature Tj. At 
the time t = a quench to a lower temperature Tp is performed. On top of this primary 
relaxation process a second one is activated by switching on an external field at the time 
t w > 0. Characterizing the response of the system under the action of the perturbation 
by the response function R(t,t'), the search for a fluctuation dissipation relation aims to 
connect R(t,t') to the relevant correlation function C(t,t') in the unperturbed relaxation 
process. Given R(t,t') and C(t,t') one can always write 

Wtf = ^rp--j$C{t,l!) (1) 

where t > t' . Without any further specification this is just a definition of the quantity 
X(t,t'), which is called the fluctuation dissipation ratio (FDR). Eq. (|J) acquires predictive 
power when independent statements are made about the FDR. Thus, if the shortest time t' 
is greater than the equilibration time t eq , i.e. if one looks into the time translation invariant 
equilibrium dynamics, the fluctuation dissipation theorem (FDT) requires X(t, t') = 1. This 
is no more true when t' < t eq . However, if appropriate conditions on t',t,t eq are satisfied it 
may turn out that X(t,t r ) depends on the time arguments only through C(t,t'). This was 
first discovered by Cugliandolo and Kurchan [p]] in the context of mean field spin glass 
models at low temperature. In that case the system does not equilibrate (t eg = oo) and 
X(t,t f ) = X(C(t,t')) in the asymptotic limit of large times. Subsequently the validity 
of this relation was verified for finite dimensional spin glass models @ and also in the 
coarsening processes of non disordered systems @,f|. As a matter of fact a classification of 
slowly relaxating systems can be made |J on the basis of the behavior of X(C). 

The relation between X and C is important for different reasons. From the point of view 
of analytical calculations it allows to close the equations of motion for R and C [|TJ. From 
the more fundamental point of view of the understanding of the off equilibrium dynamics 
it can be related to the effective temperature of different dynamical modes |J and under 
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certain hypothesis it provides a connection between the relaxation regime and the structure 
of equilibrium states J7[. 

In this paper we analyse the relaxation to equilibrium in the one dimensional kinetic 
Ising model quenched from the initial temperature Tj to the lower final temperature Tp. 
The correlation length then grows from some initial value £/, which we assume 0(1), to the 
final value 

^ = -[logtanh(J/T F )]- 1 (2) 

where J > is the nearest neighbors ferromagnetic interaction. The equilibration time is 
defined by 

By lowering the temperature of the quench t eq can be tuned at will with \im.T F ^ot eq = oo 
allowing for the investigation of the slow relaxation coming from the high temperature side. 
We compute the response function R(t, t') after switching on a random external field at the 
time t w after the quench. We are then able to analyse in detail the changeover from the 
equilibrium to the off equilibrium regime by monitoring the change in the FDR (or in the 
integrated response) as t w is varied from t w > t eq to t w < t eq . When t w > t eq dynamics is 
time translation invariant and the usual FDT holds. When the region t w < t eq is entered 
the deviation from FDT occurs. However, if the difference between t w and t eq is sufficiently 
large, there is a range of values between t w and t eq where C(t, t') scales. Whithin this range 
the FDR depends only on C, while outside there remains an explicit dependence on t w . In 
other words with a finite but large t eq the off equilibrium dynamics follows the pattern of 
interrupted aging and we find that the FDR depends only on C as long as aging holds. The 
case of the zero temperature quench is the limiting case (t eq = oo) where aging occurs for 
arbitrarily large times yielding X(t,t') = X(C(t,t')) at all times. 
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II. UNPERTURBED CORRELATION FUNCTIONS 



In the following we consider a one dimensional ferromagnetic Ising model with nearest 
neighbor interaction 

JV 

H[<r} = -Jj2wi+i (4) 

i=i 

evolving in time with Glauber single spin flip dynamics 

^P([a],t) = Y.M-a t )P([R t a],t) - w^PiHt)} (5) 

where P([a],t) is the probability of realization of the configuration [a] at the time t, [Ri<j] 
is the configuration with the i-th spin reversed 



1 ~~ ^iiPi-X + 



(6) 



is the transition rate from [a] to [Ri<j] and 7 = tanh( 



2J- 



F 



Given an initial probability distribution P([a],t = 0), with the choice @ for the 
transition rate the solution of (|5]) for large time reaches the equilibrium Gibbs state 
-P e? [c"] = \ exp(— 7=-T~t[cr])- The dynamics in the one dimensional case has been solved 
by Glauber ||. Let us summarise those properties of the time dependent correlation func- 
tions which will be needed in the following. Assuming that the initial state -P([cr], t = 0) is 
symmetrical the probability distribution P([a],t) remains symmetrical throughout yielding 
< <Ji(t) >= for all time. The equal time and the two time correlation functions then are 
defined by 

D ij (t) = ^a i a j P([a],t) (7) 

M 

C rj (t,t') = £ a lC r>P({a>},t')P(lalt' \ [a],t) (8) 

MM 

with t > t' and Cij(t,t) = Dij(t). P([a'],t' | [&],€) is the conditional probability to find the 
system in the configuration [a] at the time t, given that it was in the configuration [a'] at 
the earlier time t'. We assume that space translation invariance holds at all times. 



From (|5|) and (0) one can show that the equal time correlation function satisfies the 
equation of motion 

jD^t) = -2£> y (f) + | + Aj+i(*) + + A+i,^)] (9) 

for i ^ j and 4 = 0, since Du(t) = 1. Similarly, from (§) taking into account that 

also the conditional probability satisfies the master equation (H) one finds 

J^.(t) = -cv(m') + 1 [<Wt,f ) + c i+1>j (t,t')) . (io) 

The single spin conditional expectation can be computed explicitely obtaining || 

EPHflH*v i = i;^ w (t-o (ii) 

where -Fj_ m (t — t') = e ~ /j-^, (7 — t')) and J„(x) are the Bessel functions of imaginary 
argument. Then, using the definitions (0) and (H), the two times and the equal times 
correlation functions are related by 

C tl (t,t') = Y,D ll (t')FUt-t') (12) 
i 

or in Fourier space one finds 

C k (t,t') = D k (t')e-^-^ (13) 

with 7^ = 1 — 7 cos A;. Given the initial condition Aj(0) Eq.(|9|) can be solved exactly ||. 
After some microscopic time to, which we assume much smaller than t eq , memory of the 
initial condition is lost and for fc < 1, ^j? > 1 one has |§ 

Dk (t>) = 2(-) 1 /a 1 _ /' d^-^^e-^V,, + jfc 2 e (-[* a +€^]f+* 2 f»)l ( M ) 

7T K -\- £p JO 

where we have expanded 7^ to lowest order in k and ^p 1 using 7 = 1/ cosh(^ 1 ). The form 
(ini) for the equal time structure factor obeys the scaling relation Dfc(t') = (t') 1 ^ 2 g(k 2 t f , t' /t eq ) 



with the limits 



, ite q ) ll2 9eq{kH eq ) for t'lt e a > 1 

D k (t') ~ { V 97 " V 97 7 9 (15) 

(t') 1/2 9sc(k 2 t') for t'/t e5 « 1. 



Inserting ([TJ]) into (|T3|) and inverting the Fourier transform, the corresponding scaling form 
for the two time correlation function is obtained Cjj(i, t') = f ( Im^/l ? with the limiting 

behaviors 



, fen (iNm, 7=H fOT *7*eo > 1 

/«c(^,y) for t'/te, < 1. 
The use of the small k approximation ([14]) is justified since for small enough Tp the structure 
factor builds up a large and narrow peak about k = which gives the main contribution 
to the integral over k. In the following it will be sufficient to consider the autocorrelation 
function (i = j). In the case of the zero temperature quench the scaling behavior at the 
bottom of ( |i6|) is obeyed for all times since t eq = oo and the explicit form of the scaling 
function is given by P,|10| 



CiAt, t ) = fscit/t ) = - arcsiii 

7T 



With Tp > and t eq < oo, Ci^(t, t') can be computed by a combination of analytical and 
numerical integration (Appendix). In Fig.l we have plotted logC^^t, t') against x = (4 — 1) 
for different values of r = t'/t eq illustrating the crossover from the scaling form (|T7|) to the 
exponential decay corrasponding to the top of (|T6|) as r grows from small to large values. 
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III. RESPONSE FUNCTION 

As stated in the Introduction, let us now assume that after the quench to Tp, at some 
time t w > 0, a site and time dependent external field hi(t) is switched on. We are interested 
in the response in the magnetization to the action of this field. More specifically, we wish 
to investigate the relation between the magnetization response and the correlation function 
in absence of the field. 

For sufficiently small external field, the response in the magnetization at site i and t > t w 
is given by linear response theory 

A < <Ti(t) >=< Oi(t) > h -< ai(t) > h=0 = £ f dt'R^it^h^t') (18) 
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where 



Rij(t,t) 



5 < <Ti(t) > h s 



(19) 



8hj(f) J h=0 

is the causal response function. The difference between the two expectation values in fll8|) 
is given by A < a { (t) >= £ [ct] °AP(H *) where AP([<r],t) = P^H*) - P([a],t) is 
the difference between the probabilities with and without the field. The time evolution 
of Ph([a],t) is given by the master equaltion (H) with the transition rate (H) replaced by 
Wh(ci) = w(<7i) + Aw(<Ji) where Aw(ai) = — tanh(^)o"jw(crj). Taking hi/Tp sufficiently 
small tanh(^) ~ ^- and following || up to first order we have 

AP([a],t) = ^EE^ f dt 'Ht') [w(a' i )P([a'],t')+ 

^ F r_/i Jt w 



w(-o§P{[Ri^l!))P{W\J \[a],t) 



(20) 



yielding 



P([a'],t'l 

Performing the sum over [<r] first and using flllD we find 



(21) 



i r - 



t'1 i 



-a')p([py],olp«(*-*' 



(22) 



and since only the term with / = j survives in the summation we have 



Ri,j(t, t) 



1 " \{Dw-i{t') + D j>j+1 (t')) F^(t - t'). 



(23) 



In order to recast this result in terms of Cij(t, if) let us differentiate (|12"D with respect to 
the time arguments 



i r C i At,t') = Y j d ^±FUt-t') 
df M Y df v ' 



(24) 



^«(*,0 = -E^(0^i-i(*-0 



(25) 



Adding and (23) the summation in ([25]) cancels the second one in (E3) yielding 



(26) 



Next, using (|9]) for the time derivative when I ^ j, adding and subtracting a similar contri- 
bution with I = j and using translational invariance we can rewrite 

2 " ]\D h3+ i{t') + ^_i(0]} *<_.,■(* - 0- (27) 

Using ( |T2"D and ( [TO]) the first sum in the right hand side is given by 2&Ci ! j(t,t'), while 
the second term coincides, up to a constant factor, with the right hand side of (|23|), since 
Djj(t') = 1. Therefore, we finally get the following expression for the response function 

1 



Ri,j (t, t) 



2Tp 



d d 
— C itj (t, t ) - —Cij(t, t ) 



(28) 



which can be rewritten in the form ([!]) 



Ri,j{t,t) = — '— —C itj {t,t) 



(29) 



with 



Xi,j(t, t') 



M 
at 



~Cij(t, t') 



(30) 



Alternatively, we can expose the deviation from FDT through an additive term 

Id ,1 

Ri,j(t,t') = ^r-Q^C itj (t,t) - ■^rB iij (t,t') 



(31) 



with 



, d d 
B id (t,t) = —C id (t,t') + —C id (t,t'). 



(32) 



dt> t>JV ' ' • at 

When time translation invariance holds we have either Xij(t,t') = 1 or Bij(t,t') = and 
the usual FDT is recovered. 
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IV. RANDOM EXTERNAL FIELD 



In the simulations [|J the external field is taken random with site independent bimodal 
distribution 



P[h] = IT, 



-5{hi -h) + U{h t + h) 



(33) 



The reason for this choice is not to bias the evolution toward the formation of predominantly 
positive or negative domains through the introduction of the external perturbation. The 
quantity of interest then is the staggered magnetization 

1 



N 



hi 
h 



(34) 



where the bar represents the average over the field configurations. From flTB]) and (^) the 
integrated response is given by 



Xu(t,t w ) = ]-M(t,t w ) = f dt'R u (t,t r ). 
n Jt w 



(35) 

Dropping the double index and inserting the form (p9|) of the response function, if the 
FDR depends on time only through C(t,t') we have 



T FX (C{t,t w )) 



C(t,t w ) 



dCX{C) 



(36) 



namely also the integrated response depends on time only through the autocorrelation func- 
tion. This occurs when FDT holds with X(C) = 1 yielding 



T FX {C{t,t w )) = [l-C{t,t w )\ 



and when the scaling form (fPTD holds. In that case from (|30| ) follows 



(37) 



x(t,t') = - 



(38) 



and inverting (IT 



l> sin" 



(fCfttf) 



t 2-sin 2 (fC(t,t'; 
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(39) 



we find 



X(C) 



1 



(40) 



2 — sin 




Inserting this into (|36|) we obtain 



Tpx(C(t,t w )) = — arctan \/2cot(— C(t, t w )) 



(41) 



We have then proceeded to compute (Appendix) TFx(t,t w ) with the values of the pa- 
rameters (t eq = 10 3 ) corresponding to the behavior of the autocorrelation function dis- 
played in Fig.l and we have plotted Tpx(t,t w ) against C(t,t w ) in Fig.2 for different values 
of t w /t eq . In order to understand the plot notice that if t eq is finite, from ( |5o] ) follows 
lim^oo Tpx{t, tw) = ^-M eg = 1 where M eq is the equilibrium value of the magnetization. 
On the other hand, with a finite t eq one has also lim^oo C(t, t w ) = 0. Therefore, when 
plotting TpX vs - C all the curves starting out at (C = 1, TpX — 0) must end up in the same 
point (C = 0,Tpx — !)• The dependence on t w /t eq enters on how the initial and the final 
point are joined. Thus, if t w /t eq > 1, FDT holds over the entire time interval (t w ,t) and 
the plot is linear according to (p7|). However, if t w /t eq < 1 then it is possible to have also 
t/t eq < 1. In that case C(t,t') obeys the scaling form (|T7|) and in the range of values of C 
where this holds, TpX follows the shape (|4~1|). This forces the plot to fall below the straight 
line of the FDT, but eventually as C decreases the plot must raise again in order to reach 
the value T F x = 1 at C = 0. Therefore, the peculiar shape of the curves displaying a change 
in concavity is a consequence of a finite equilibration time. The final upword bending of 
the curves corresponds to interrupted aging and that is where the curves do depend on t w . 
Furthermore, the range of values where the plot follows the shape (^Tj) is larger the smaller 
is the value of t w /t eq . In the limiting case t eq = oo aging holds for all time and the plot 
obeys d4l|) over the entire range of C values. 
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V. CONCLUSIONS 



The relaxation dynamics of the one dimensional Ising model allows to analyse in detail the 
transition from the off equilibrium to the equilibrium regime. In particular, we have obtained 
the crossover in the FDR from the nontrivial form X(C) given by ([5]) to X(t,t') = 1 as 
a manifestation at the level of the response function of the crossover in the underlying 
correlation function from aging to time translation invariance. 

A comment should be made about the shape of X(C). In the case of the zero temperature 
quench ( ffij ) holds for all time. On the other hand, the zero temperature quench is a phase 
ordering process eventually leading to the coexistence of ordered phases as in the quench 
below the critical point of a system with a finite critical temperature. In the latter case 
X(C) displays ||,[5] a qualitatively different behavior decreasing from 1 and flattening to 
zero, while in our case X(C) decreases from 1 toward 1/2 as C goes to zero. Although we 
do not have a complete understanding of the origin of this discrepancy, we believe this to 
be related to the absence in the one dimensional case of the asymmetry term in the relation 
between R(t, t') and C(t, t'). In the context of Langevin dynamics one can derive [11 1 in full 
generality 

^'M^-^h^-m^ (42) 

where A(t, t') is the asymmetry term which vanishes for linear dynamics. From fli2"|) the 
FDR takes the following general form 



C(t,f) 



1 A(t,t') 



(43) 



2£C(t,tf) 



and if we assume scaling C(t,t') = f{t/t') the square brackets contribution is given by ( |38|) 
independently from the form of f(x). Therefore, if we accept that A is a function of C when 
scaling holds, in order to have limc^o X(C) < 1/2 the asymmetry term must necessarily be 
nonzero. Now, from Eq.(p8|) follows that in the one dimensional Ising model with Glauber 



dynamics the asymmetry is absent. Indeed, in this case as Eq.s (pT) and (|T0D show, dynamics 
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is linear. Another example of linear dynamics leading to limc^o X(C) = 1/2 is the massless 
gaussian model []TT] . Conversely, if one considers the Ising model with higher dimensionality 
and a finite critical temperature, the equations of motion for the two point correlation 
functions are coupled to higher correlation functions producing a nonlinearity which in turn 
is expected to produce a nonvanishing asymmetry in the off equilibrium regime. 



VI. APPENDIX 

In order to carry out the computation of Tpx(t,t w ) we start from the sum of the time 
derivatives of Ck(t,t') obtained from (O) 



Using (14) for Dh(t') and carrying out integrations by parts we find 
dtf V 7T Jo Jy 



Inserting this result into (|44]) and integrating over k we have 



d d\ . dkdD k {t') 



B(M ' ) = (aF + s) c(i ' i ' ) = / 



2tt dtf 



e -\{^+^- 2 ){t-t') 



I 7T JO 



1 _ \2(t - f) ^u-2 {t+tl) 



ir(t + V) V t' 

Next, Fourier transforming the equation of motion (^j) one obtains 



with 



e 2 < =f L 



(45) 



TT\t'(t-t') 7T Jo V !/(* + *' - 2yt') 3 

e -|^ 2 (*+*'). ( 46 ) 



-££> fc (*) = -2 7fc D fc (t)+r(t) (47) 



r(t) = — =^ + ^r/( > /^ a t) (48) 

V 7TC 



where £Jr/ is the error function. Inserting (|47D in the right hand side of fl44]) and using 



|C fe (t,t') = -7*A(M') one finds 
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'd__cT 
dt 1 dt 



C k (t,t') =r{t')e-^ k (49) 



which after integration over k gives 

at; *J2t'(t-t>) U Vt^v 



Inserting this result in fl28|) with z = j the integrated response Tpx(t, t w ) is obtained carrying 



out numerically the integration in ( |3~5"D 



Similarly, the autocorrelation function is obtained by taking the difference of ( |4~6|) and 
(|50|) and carrying out numerically the time integration in 

C(t,t w )=l + J* ds^C(s,t w ). (51) 
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Figure Captions 

Fig.l - Plot of the autocorrelation function for different values of r = t'/t eq and t eq = 10 3 . 
The continuous line is the plot of the anlytical solution ( [FTP coresponding to Tp = 0. 

Fig. 2 - Plot of the integrated response for different values of t w /t eq and t eq = 10 3 . The 
continuous line is the plot of the anlytical solution (fl"H ) coresponding to T F = 0. 
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